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Neutron penumbral imaging technique has been successfully used as the diagnosis method in Inertial Confined 
Fusion. To help the design of the imaging systems in the future in CHINA. We construct the Monte carlo 
imaging system by Geant4. Use the point spread function from the simulation and decode algorithm (Lucy- 
Rechardson algorithm) we got the recovery image. 



Inertial confinement fusion (ICF) has attracted much 
attention from all over the world. The National Ignition 
Facility (NIF) programs in America and Laser MegaJoule 
(LMJ) in France aim at reaching the inertial confine- 
ment fusion of a deuterium-tritium (DT) filled cryogenic 
target. China's ICF program are also in progress, this 
paper's aim is just study the neutron diagnosis system 
which will be used in the future ICF program in China. 

The process of ICF can be divided into four steps as 
Fig 1 : First step generate radiation to drive implosion. 
There are two approaches to drive implosion, one is laser 
direct driven, another is laser energy injected inside the 
cavity generates X-rays that heat the DT fuel. Under 
X-ray or laser compression, the target density as well as 
the mean ions' temperature increases. Fusion reactions 
will occur within a hot spot generating 14 MeV neutrons 
and alpha particles^ If the density is sufficient to trap 
the energetic alphas or other radiation, these energy will 
be deposited on the surrounding cold DT medium. Then 
the chained thermonuclear reactions will lead to a com- 
bustion in the compressed target. For a well compressed 
target, combustion will last long enough to generate more 
energy than the injected energy. If ever, drive symmetry 
is poor or pulse shaping is imprecise, burning will be low. 
So we can use the neutrons escaped from the compressed 
cores of target to study the compressed process to make 
sure the drive be symmetrical. 

There are several neutron imaging techniques for in- 
ertial confinement fusion (ICF) experiment to detect the 
size, shape, and uniformity of the compressed cores. The 
pinhole imaging is the easiest way, but the low neutron 
yield has an effect on the sensitivity. In NIF experiment 
a pinhole array is designed to increase the sensitivity^^. 

The another technique is the neutron penumbral imag- 
ing (NPI), which have been widely studied by LMJ^^^ 
and other groups^^^^. Compared with other neutron 
imaging techniques, the sensitivity of neutron penumbral 
imaging is higher. Annular imaging is nearly the same 



^Wangphysics@126.com 
^)tangyq@126.com 




FIG. 1. ICF physics process (Radiation compressed, Implo- 
sion, Hot spot creation, Combustion) 



as penumbral imaging, the only difference of the imaging 
process is aperture. 

In this paper we just based on the simulation model 
in the reference^^^^ to design a software package NG- 
PINEX which can be used to simulate the neutron imag- 
ing process and recovery the coded image. The imag- 
ing process simulation have been described in Geant4, 
and then the recover algorithm we used Lucy-Richardson 
method^"^'^^, which have been successfully applied on de- 
coding x-ray ring code image^^. The whole program is 
written in C++, which can be run under the linux system 
with multi-core server. 

In the program we simulated neutron penumbral imag- 
ing and Annular imaging, the pinhole imaging block is 
also included in the program. Based on the NIF and 
LMJ experiment we design an imaging system aimming 
to achieve 20fim resolution. At last we use the Lucy- 
Richardson method successfully got the recovered image 
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FIG. 2. Penumbral imaging 



with 20/im resolution 



I. THE IMAGING SYSTEM SIMULATION 

The structure of neutron penumbral imaging system 
is shown in Fig. 2. The source , aperture and detector 
composed the whole imaging system. The fast neutrons 
emitting from a burning target are scattered by the aper- 
ture materials, and form an penumbral image on the de- 
tector arrays. In our simulation the coded aperture is 
placed at 55mm from the source, the detector is placed 
at 7m from the aperture. The aperture is designed as 
a 50mm{h) thick 10mm diameter(tungsten) cylinder to 
achieve an effective absorption of the lAMeV neutrons 
in the opaque parts of the aperture. The pattern of the 
aperture is a disk (penumbra) which consists in a biconi- 
cal hole, the diameters (I^i, D2-, Ds) of the biconical hole 
are 0.3mm, 0.38mm, 0.535mm. An enlarged coded im- 
age is projected on the neutron detector which is placed 
on the target bay floor, at 7895mm from the aperture. 
The image obtained from the detector is a convolution of 
the source spatial distribution with the aperture transfer 
function. 

In the Fig. 2 only the red annular of the coded im- 
age contains the effective information of the source shape 
changing. To achieve the high SNR we can place a plug 
inside the hole of an aperture of the same dimensions 
as the penumbra(Fig.2) to define a annular, which will 
be study in another paper. This pattern can shield a 
lot of background but also face the alignment tolerance 
problem and not enough neutron can be detected. 

The neutron detector is composed of 255 x 255 array 
of scintillating fibers Fig. 3. Each scintillating fiber can 
be realized as a sub-detector, the structure was designed 
as the left figure of Fig. 3. The core part is a tubes of 
460/im diameter and the material is polystyrene BCF- 
10 (molecular formula C6H5CH=CH2, img/cm^). 14 
MeV neutrons mainly interact by elastic scattering on 
hydrogen nuclei inside the fibers. As they lose kinetic 
energy, the recoil protons will produce light, a part of 
the light is guided through these fibers. Outside the core 
part is a layer of the clad material, which is used to avoid 



Clad 



C ore 




FIG. 3. Fiber array 



the scintillating light interaction in the array of scintil- 
lating fibers. The material of the envelope is Polyvinyl 
alcohol (molecular formula CH2=CH0H, 1.26^/cm^) and 
the thickness is 20/im. An optical relay then casts the 
scintillating light from the end of the array to an image 
intensifier tube (IIT) which gates the 14 MeV neutrons 
at their arrival time on the detector. The image is then 
reduced with a fibered optical taper onto a CCD. Overall 
spatial resolution of the system for both the annular and 
the penumbral apertures is the same, which is determined 
by the point spread functions (PSF) of the coded aper- 
ture and detector resolution {/S.s detector) scaled down the 
source plane as equation 1. In the paraxial approxima- 
tion, the broadening of the apertures point spread func- 
tion is identical, determined by the FOV and the aperture 
to target distance (i^o)- 

\ln(2) X FOV 7^ 

where ji is the attenuation of the neutrons in tungsten 
and Li is the detector to the target distance. With 
this design, the overall spatial resolution of the system 
is about 20 /im in source plane. 

In the GEANT4 simulation we used the physics list 
QGSP_BERT_HP which contains the high precision neu- 
tron package (NeutronHP) to simulate the process of 
neutrons' transportation with the energy below 20 MeV 
down to thermal energies. The physical process neutron 
capture, fission, elastic and inelastic scattering (including 
absorption) are treated by referring to the ENDF-B VI 
cross-section data. Besides, in the ref^^ has demonstrated 
the validity of GEANT4 calculations in neutron gener- 
ation and transportation which is the same as MCNP. 
Moreover, object-oriented programming and highly so- 
phisticated processing of both electromagnetic interac- 
tions and ionization processes enhance GEANT4 capa- 
bilities for the study of mixed fields and complicated ge- 
ometries. 

In the ICF experiment the neutrons generated by hot 
spot spread in whole space. The idealized aperture only 
allow the neutrons in the solid angle as Fig. 2 to deposit 
the energy on the detector. But actually because the neu- 
trons have the high ability to penetrate the start part of 
the aperture, and also the neutrons will scatter with the 
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FIG. 4. The neutron radiated region in the simulation 



wall of the aperture, these neutrons which hit the de- 
tector will become the background of the imaging. To 
simulate the complete imaging process and also save the 
calculation time we ask the program only generate ran- 
dom distribution neutron in the cone with Ri radius as 
Fig. 4. As equation 2 Ri is determined by Rq which is 
the width of the detector project on the front face of the 
aperture. 

With this design there are about one half of the events 
be scattered by the penumbral aperture. To get the high 
statistic sample we simulated 1 x 10^^ effective events of 
penumbral imaging. 

r Do/2, Ro > Do/2 
Ri = <^ 1.3i?o, Ds/2 <Ro< Do/2 (2) 
[ Ds/2, Ro < Ds/2 

II. RECONSTRUCTION OF THE CODED IMAGE 

At present, there are a lot of reconstruction methods 
have been developed. Among them the simplest and 
quickest method is Wiener filter method, which is a linear 
reconstruction method. As we cannot get prior knowl- 
edge of the signal-noise ratio, wiener filter method cannot 
get good results. Therefore we can use this method as the 
reference or use the result as the input of some nonlin- 
ear reconstruction method. In these years there are some 
nonlinear reconstruct methods have been introduced to 
reconstruct penumbral image such as genetic algorithm 
and heuristic algorithm proposed by Chen et al.^^, and 
the classical molecular dynamics reconstruction method 
proposed by Liu et al.^^'^^. 

In the experiment L. Disdier et al. had success- 
fully used the filtered autocorrelation technique^^ to get 
20/im resolution reconstruction image from the penum- 
bral imaging system with 22.3/im theoritical resolution. 
Filtered autocorrelation technique is different from other 
method which don't need to rely on the simulation PSF 
in the reconstruction. 

Another nonlinear method Lucy-Rechardson 
algorithm^^'^^ has been successfully used on X-ray 
ring coded imaging in experiments^. And this algorithm 
has been most widely used for restoring images with 
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FIG. 5. Distribution of point spread function(PSF). 




FIG. 6. 'F' source imaging for penumbral aperture (left), un- 
folded result (right). 



noise from counting statistics and the data approximate 
Poisson statistics. In this program we first try to use the 
Lucy-Rechardson method to get reconstructed image. 
The process of the imaging can be described as a spatial 
intensity distribution g{x,y) convolved by a neutron 
source spatial intensity distribution /(x, y) with a point 
spread function (PSF) h{x,y). Add the noise n{x,y) the 
detector image can be expressed as 

g{x, y) = f{x, y) (g) h{x, y) + n(x, y) (3) 

Lucy-Rechardson algorithm is an iterative method. 
Based on the reference^^'^^, the decode formulation can 
be expressed as 4. Where PSF(h) and image on detec- 
tor (g) come from the simulation, the reconstructed im- 
age (f) can be extracted after k times iterate. 

fk+i = fk{hQj-^) (4) 

The PSF for penumbral and annular imaging are pre- 
sented in Fig. 5 respectively. In Fig. 6 we presented the 
raw images with a 'F' shape neutron source and the un- 
folded results with Lucy-Rechardson method. The width 
of the 'F' shape source is 20/im. 

From the Fig. 6 we can see with this design the imag- 
ing system can achieve the aim to reconstruct the neutron 
source. With Lucy-Rechardson method the unfolded im- 
age of Penumbral aperture is clear. To get the real res- 
olution of the system we simulate two point source with 
distance 20/im imaging on the detector. In the Fig. 7 
present the raw image and unfolded image. For the an- 
nular imaging reconstruction and detail of the unfolding 
method we will present in further study. 
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FIG. 7. Two point source imaging for penumbral aper- 
ture(left), unfolded result (right). 



III. CONCLUSION 

We have successfully construct one software NG- 
PINEX which integrated the imaging system simulation 
and the reconstruction method. The simulation and re- 
construction result show we can use this imaging system 
with these parameters to get 20/im resolution image of 
DT implore. To attend the aim of 5/im resolution in the 
future we can use this software to optimize the system 
parameter. 
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